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ABSTRACT 

The whole set of INTEGRAL observations of type la supernova SN2014J, covering the period 19-162 
days after the explosion has being analyzed. For spectral fitting the data are split into early and late 
periods covering days 19-35 and 50-162, respectively, optimized for 56 Ni and 56 Co lines. As expected 
for the early period much of the gamma-ray signal is confined to energies below ^200 keV, while for 
the late period it is most strong above 400 keV. In particular, in the late period 56 Co lines at 847 and 
1248 keV are detected at 4.7 and 4.3 a respectively. The lightcurves in several representative energy 
bands are calculated for the entire period. The resulting spectra and lightcurves are compared with 
a subset of models. We confirm our previous finding that the gamma-ray data are broadly consistent 
with the expectations for canonical ID models, such as delayed detonation or deflagration models 
for a near-Chandrasekhar mass WD. Late optical spectra (day 136 after the explosion) show rather 
symmetric Co and Fe lines profiles, suggesting that unless the viewing angle is special, the distribution 
of radioactive elements is symmetric in the ejecta. 


1. INTRODUCTION 


A Type la supernova is believed to be a therm onuclear 
explosion of a carbon-oxygen (CO) white dwarf (|Hoyle fe 
Fowler| 1960) in a binary system, (see, e.g.,|Hillebraudt fe 


IN iemeyer| 21)00 [ |Imshennik fe Dunina-Barkovskaya 


TEGRAL has been reported in |Churazov et ah (|2014b|) 
(days ~5 0-100 since explos ion), IDiehl et al.|(|2014|) (days 
^16-1 9 ), llser n et al.|(120151) (days ~ 16-35), see aisolDiehll 
et al. (2015 k Despite of the proximity, SN2014J in 


for a review). Most popular scenarios of the explo- 


Chandrasekhar limit (e.g., 

Whelan & Iben 1973), (ii) a 

merger/collision of two WDs (e.g., 

Iben & Tutukov|1984 

Webbink 1984 

Kushnir et al. 2013 

), (iii) an initial explo- 


gamma-rays is an extremely faint source and the ex¬ 
pected signal is below 1% of the background. This makes 
the results sensitive to the adopted procedure of the 
background handling by different groups and lead to ten¬ 
sion between some results. Here we have combined all 
INTEGRAL data and uniformly proce s s them using the 
same procedure as in Churazov et al. (2014b|). The re- 


triggers subsequent explosion of the bulk of the material suiting spectra and light-curves are compared with the 

v f-\ r \ inti Ann rvr Knm A Ft ryx a I A vn Anl A n 


(e.g., Nomoto & Sugimoto 1977 Hoeflich & Khokhlov 


1996). In all scenarios a thermonuclear runaway converts 


predictions of basic type la models. 


substantial fraction of CO mass into iron-group elements 
and the released energy powers the explosion itself. The 
optical light of the supernova is in turn powered by the 
decay of radioactive elements, synthesized during explo¬ 
sion. For the first year since the explosion the decay chain 
of 56 Ni—> 56 Co—K 56 Fe is of prime importance. As long as 
the expanding ejecta are optically thick for gamma-rays 
the bulk of the decay energy is thermalized and is re¬ 
emitted in the UV, optical and IR band. After several 
tens of days the ejecta become optically thin for gamma- 
rays making SNIa a powerful source of gamma photons. 

Here we report the results of INTEGRAL observations 
of SN2014J covering a period from ~16 to ^162 days 
since the explosion. 

The analysis of the SN2014J data obtained by IN- 


plosions (e 

R-i 

Moll et al. 

201 


Seitenzahl et al. 

2013 

Fink et al. 

2014 


4| lead to a complicated distribution o 


angle dependence in the predicted gamma-ray flux. How¬ 
ever, the overall significance of the SN2014J detection 
in gamma-rays by INTEGRAL (see ^3] and f|5| corre¬ 
sponds to ~ 10 s.t.d. This precludes a very detail model- 
independent analysis. We therefore took a conservative 
approach of comparing the data with a subset of popu¬ 
lar ID SNIa models (see Q, some of which were used 


Milne et al. (2004) for assesment of SNIa gamma-ray 


codes. While these models do not describe the full com¬ 
plexity of SNIa ejecta, they can serve as useful indicators 
of the most basic characteristics of the explosion, includ¬ 
ing the total mass of radioactive nickel, total mass of 
the ejecta and the expansion velocity. We also verify 
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(^5.3) if adding an extra component, corresponding to a detection of gamma-rays from SN2014J possible. 


transparent clump of radioactive Ni, on top of the best¬ 
fitting ID model, significantly improves the fit. In £j 6 ] 
we make several basic consistency checks of gamma-ray 
and optical data, using optical observations taken quasi- 
simultaneously with INTEGRAL observations. Section 
[7] provides the summary of our results. 

2. SN2014J IN M82 


SN2014J in M82 was discovered ( 

Fossev et al. 

20141 

on Jan. 21, 2014. The reconstructed 

(Zheng et al. 

2014| 


UT with the uncertainty of order ±0.3 days. At the 
distance of M82 (~ 3.5 Mpc), this is the nearest SN 
la in several decades. The proximity of the SN2014J 
triggered many follow-up o bser vations, including those 
by INTEGRAL pCuulke rspMI ). 


The SN is located ~1 kpc trom the M82 nucleus and 
has a strong (Ay ~ 2 ) and complicated absorption in the 


UV-o ptical band (e.g., Goobar et al.||2014 ; Marion et al. 


2015| lAmanullah et al. |2014[ |Foiey et ai.||2014t |Welty et 

■-^Q 15 ||g — 


al. ||20141 |Patat et ai.||2015| |Brown et al.|2015| IKawabata 

et al.||2014[). ~ 


From the light curves and spectra SN2014J a ppears to 
be a “normal” SNIa with no large mixing (e.g., Marion et 


al.|2015 Ashall et al .|2 

detonation models. De 


2014), consistent wit h the delayed- 


dctonation mo dels. Detec tion of stable Ni (Friesen et al. 
2014 Telesco et al.||2015 ) in IR suggests high density o f 
the burning material (see, e.g., Shigeyama et al.| 1992), 
characteristic for near-Chandrasekhar WD. 

Search in X-ray, radio and optical bands (including 
pre-supernova observations of M82) didn’t reveal any ev¬ 
idence for accretion onto the WD before the explosion, 
any candidate for a companion star, or compelling evi¬ 
dence for a large amount of circumbina ry materia l, im¬ 
plicitly supporting the DD scenario ‘(Kclly~et al.|[2014] 


Nielsen et al.||2014l |Margutti et al.|2014[ |Perez- Torres et| 

al.||2014p , although some SD scenarios are not excluded. 

In gamma-rays the first detection of SN2014J in 56 Co 
lines was reported abou t 50 days since the explosion 
(Churazov et al. 2014a). The gamma-ray signal from 
S1N2014J was also report ed in the earlier phase ~16-35 


days after the explosion (Isern et al. 2014 Diehl et al. 
20141. 

Throughout the paper we adopt the distance to M82 
(and to SN2014J ) of 3.5 Mpc. The recent analysis by 


3. INTEGRAL OBSERVATIONS AND BASIC DATA 
ANALYSIS 

INTEGRAL is an ESA scientific mission dedicated to 
fine spectroscopy and imaging of celestial 7 -ray sources 
in th e energy range 15keV to lOMeV (Winkler et al. 
20031. 


The INTEGRAL data used here were accumulated 
during revolutions 1380-1386, 1391-1407 and 1419-14210 
corresponding to the period ~16-162 days after the ex¬ 
plosion. 

In the analysis we fol l ow the proc e dures described in 
Churazov et al. (2014b); Isern et al. (2013) and use the 
data of two instruments SPI and 1SGK1/1BIS on board 
INTEGRAL . 

3.1. SPI 

SPI is a coded mask germanium spectrometer on board 
INTEGRAL . The instrument consists of 19 individual 
Ge detectors, has a field of view of ~30° (at zero re¬ 
sponse), an effective area ~ 70 cm 2 at 0.5 MeV and en¬ 
ergy resolu tion of ~2 keV (|Vedrenne et al.|2003| Roques 
et al.|[2003 ). Effective angular resolution of 8 P 1 is 


72 ° 


During SN2014J observations 15 out of 19 detectors were 
operating, resulting in slightly reduced sensitivity and 
imaging capabilities compared to initial configuration. 

Periods of very high and variable background due to so¬ 
lar flares and passage through radiation belts were omit¬ 
ted from the analysis. In particular, based on the SPI 
anti-coincidence system count-rates, the revolutions 1389 
and 1390 were completely excluded, as well as parts of 
revolutions 1405, 1406, 1419, 1423 and 1426. The data 
analysis follows the scheme implemented for the anal¬ 
ysis of the Galactic Center pos itron annihilation emis¬ 
sion (IChurazov et al.|2005l 
al. 2001 


2011). We used only “single” 
events (|Vedrcnnc et al. 2003) and for each detector, a lin¬ 


ear relation between the energy and the channel number 
was assumed and calibrated (separately for each orbit), 
using the observed energies of background lines at 198, 
438, 584, 882, 1764, 1779, 2223 and 2754 keV. 

The flux of the supernova S(E) at energy E and 
the background rates in individual detectors Bi(E,t ) 
were derived from a simple model of the observed rates 
Di(E,t ) in individual SPI detectors, where i is the de¬ 


Foley et al. (|2014|) suggests the distance of 3.27±0.2 Mpc. 

typical exposure of 2000 s: 


This estimate is formally consistent with the U ~ 3.53 ± 
0.26 Mpc from Karachentsev & Kashibadze|2006 and our 

« S(E) x Ri(E, t) + Bi(E,t). 

( 1 ) 


ear in mind 

that all fluxes and normalizations of best-fitting models 
can be overestimated (underestimated) by as much as 
~ 20 %. 

The only other supernova sufficiently bright to allow 
for detailed study in gamma-rays from 56 Ni and 56 Co de¬ 
cay is the Type II SN1987A in Large Magellanic Cloud. 
In SN1987A the down-scattered hard X-ray c ontin uum 
was first se e n half a year after t he explosion flSunyae v 
et al. 1987 |Dotani et al. 1987 Sunyaev et ah " 1990), 
while 7 -ray lines of 0b Co were detected several months 


later (Matz et al. 
SN2014J is more t 


1988 Teegarden et al. 1989). While 
lan 60 times further away from us than 


SN1987A, the larger amount of radioactive 56 Ni and less 
massive/opaque ejecta in type la supernovae made the 


Here Ri(E,t) is the effective area for the i- th detec¬ 
tor, as seen from the source position in a given obser¬ 
vation. The background rate is assumed to be linearly 
proportional to the Ge detectors’ saturated event rate 
Gsatit) above 8 MeV, averaged over all detectors, i.e. 
Bi(E,t) = Pi(E)G Sa t(t) + Ci(E), where Ci(E) does not 
depend on time. The coefficients S(E),/3i(E) and Ci{E) 
are free parameters of the model and are obtained by 
minimizing y 2 for the entire data set. Even though the 
number of counts in individual exposures is low, it is still 
possible to use a plain 7 2 approach as long as the errors 
are estimated using the mean count rate and the total 

1 http://www.cosmos.esa.int/web/integral/schedule- 

information 
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num ber of coun ts in the entire data set is large (Chura- 
zov et al.l 1996). The linear nature of the mode' 

(hi 


allows 


for straightforward estimation of statistical errors. 

Despite its proximity, SN2014J is still an extremely 
faint source in 7 -rays. FigjI] shows the comparison of the 
quiescent SPI background, scaled down by a factor of 10 3 
with a sample of representative models. Two models la¬ 
beled “20d uniform” and “16-35d W7” show the models 
for the early period of SN2014J observations. The for¬ 
mer model is based on the best-fitting 3PAR model to the 
SN spectra recorded betw een 50-100 days after explosion 
(Churazov et al. 2014b), recalculated for day 20. The 
model assumes uniform mixing of all elements, including 
the radioactive 56 Ni, across the ejecta. This model at 
day 20 produces prominent 56 Ni lines near 158 keV and 
812 keV. The latter model (W7, see 0 averaged over pe¬ 
riod 16-35 days does not include mixing and it produces 
much fainter lines. Finally the “50-162d W7” model cor¬ 
responds to later observations. The most prominent fea¬ 
tures of this model are the 56 Co lines at 847 and 1238 
keV. Among all these features the 56 Co line at 1238 keV 
is located in the least complicated portion of the back¬ 
ground spectrum. 

The spectral redistribution matrix accounts for 
the instrumental line broadening estimated from the 
data, accumulated during SN2014J observations. We 
parametrize the energy resolution as a Gaussian with the 
energy dependent width 


» 0.94 (E line / 500f 115 keV. 


( 2 ) 


Compared to our previous analysis we amended the 
spectral redistribution matrix of SPI by including low 
energy tails, associated with the interactions (Compton 
scattering) of incoming photons inside the detector and 
in the surrounding material. These photons are still reg¬ 
istered as single events in the SPI data, but their energies 
are lower than the true incident energy. We used the re¬ 
sults of Monte-Carlo simulatio n of SPI energy/imaging 
response (Sturner et al. 2003) and folded-in our proce¬ 


dure of spectrum reconstruction described above. For 
steep spectra the account for low energy tail results in 
a modest ~ 10 % change in the spectrum normalization, 
while for the very hard SN2014J spectrum it produces 
a low energy tail which provides large contribution to 
the continuum, while fluxes of narrow lines remain unaf¬ 
fected (Fig{2j). With this response matrix the Crab Neb¬ 
ula spectrum, observed by INTEGRAL made between 
Feb 21 and 23, 201 4, is we ll d escr i bed by a b roken power 
law obtained by Jourdain & Roques (2009) for earlier 
Crab Nebula observations with INTEGRAL . 

In our analysis we usually ignore the part of the spec¬ 
trum at energies higher than 1350 keV, since in the en¬ 
ergy range between 1400 and 1700 keV the instrument 
suffers from the enhanced detector electronic noise, while 
at even higher energies only we aker lines from 56 Co decay 
are expected (see Table [3| in | ]_1.2[ ). The convolution of 
the fiducial SNIa mod el (see jQ with the simulated SPI 
response (Sturner et al.||2003) confirmed that the contri¬ 
bution of high energy lines is negligible below 1350 keV, 
at least for “single” events considered here. 

The inspection of Fig[l]shows that there is no chance to 
detect continuum in the SPI data for any of our fiducial 
models. E.g., for a 100 keV wide energy bin between 600 


and 700 keV the expected S/N after 4 Msec observation 
between days 50 and 162 is ~ 0.5 ct. In the real data 
no evidence for significant continuum abo ve 5 00 keV was 
found in the time-averaged spectra (see j j5.1| below). As 
Fig- 0 the off-diagonal tail of the 847 and 1238 keV lines 
dominates over intrinsic SN continuum (see Fig. [ 2 ]), while 
the line shapes and fluxes are not affected. 

In general, we consider the inclusion of the off-diagonal 
term in the response as an improvement compared to a 
pure diagonal response. We used this improved r espo nse 
throughout the paper and at the same time in §5.1| we 
consider several data sets, which include or exclude the 
SPI data below ~400 keV. Inclusion of the low energy 
(< 400 keV) data boosts the S/N, while the exclusion of 
these data (dominated by off-diagonal continuum) makes 
spectral fits less prone to possible uncertainties in the off- 
diagonal term calibration. 

To verify the whole SPI pipeline, we have done an in¬ 
dependent analysis of the same data using the tools and 


procedures orig 

inallv developed and tuned for SN2011fe 

(see Isern et al. 

2013 

). This analysis includes energy cal- 
modeling and the background and 

ibration, background 


source fluxes fitting. Verification of these steps is impor¬ 
tant since the source (SN2014J) is very faint and even 
subtle changes in the calibration might result in signifi¬ 
cant changes in the source spectrum. The fluxes in the 
835-870 keV band were derived using these two inde¬ 
pendent pipelines for every revolution during SN2014J 
observations. Comparing fluxes point by point, we have 
found very good agreement, with the scatter well within 
statistical errors. The signal from SN2014J is seen in 
both pipelines. No systematic trends of deviations with 
the variations of the flux level are found. We have con¬ 
cluded that the results are fully consistent, within the 
assumptions made on the background parameterization. 


3.2. ISGRI/IBIS 

The p rimary imaging instru ment inboard INTEGRAL 
is IBIS (Ubertini et al. 2003) - a coded-mask aperture 
telescope w ith the Cd'Te-based detector ISGRI (Lebrun 
et al.|200 3). It has higher sensitivity to continuum emis¬ 
sion than SPI in the 20-300 keV rang(0 and has a spa¬ 
tial resolution ~ 12'. We note here, that neither ISGRI, 
nor SPI can distinguish the emission of SN2014J from 
the emission of any other source in M82. In particu¬ 
lar, M82 hosts two ultra-luminous and variable sources 


(e.g. |Sazonov et al.||2014| |Bachetti et al. 1 120 1 T[) which 
contribute to the flux below ~ 50 keV. 1SCR1 however 
can easily differentiate between M82 and M81, which are 
separated by ~ 30'. The energy resolution of ISGRI is 
~10% at 100 keV. The ISGRI energy calibration uses 
the procedure implemented in OSA 10.039. The images 
in broad energy bands were reconstructed using a stan¬ 
dard mask/detector cross-correlation procedure, tuned 
to produce zero signal on the sky if the count rate across 
the detector matches the pattern expected from pure 
background, which was derived from the same dataset 
by stacking detector images. The noise in the resulting 
images is fully consistent with the expected level, de¬ 
termined by photon counting statistics. The fluxes in 
broad bands were calibrated using the Crab Nebula ob¬ 
servations with INTEGRAL made between Feb 21 and 


2 http://www.cosmos.esa.int/web/integral/aol3 
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Fig. 1.— SPI quiescent background in compariso n wit h the repr esenta tive model spectra. SPI background is multiplie d by a factor 10“ 3 . 
Green and blue lines correspond to the W7 mo del jNomoto et al. 1984 b averaged over the early and late periods (see 5j3.3|) , respectively. 
Red line shows the model the 3PAR model from |Churazov et al 1 fit) 141ij for day 20 since the explosion. In this model all elements, including 
radioactive isotopes, are mixed uniformly over entire ejecta. The robust prediction of all plausible models is the presence of two 56 Co line 
at 847 and 1238 keV during the late phase. Vertical lines show two energy bands used for making images. The “cleanest” SPI background 
is near the 1238 keV line, where no strong instrumental lines are present. 


23. The Joui'dain & Roques (2009) model was assumed 
as a reference. 


3.3. Lightcurves, Spectra and Images 

The lightcurves in several energy bands were gener¬ 
ated using IGSRI and SPI data. The time bins (^3 days 
each) correspond to individual revolutions of the satel¬ 
lite. Finer time bins are not practical given that the 
source is very faint. The lightcurves are shown in Figs]3l 
[| together with a set of representative models (see (j4||. 
For the broad 100-200 keV band the conversion of the IS- 
GRI flux using Crab spectrum as a reference is not very 
accurate because of the difference in the shape of the 
incident spectra. The conversion factor has been recal¬ 
culated using several representative SN models, resulting 
in a modest ~13% correction factor, applied to the fluxes 
shown in Fig. [3] 

In principle, the spectra can be extracted for any in¬ 
terval covered by the observations, e.g., for individual 
revolutions, as is done above for the lightcurves in sev¬ 
eral broad bands. For comparison of the observed and 
predicted spectra we decided to split the data into two 
intervals covering 16-35 and 50-162 days after the explo¬ 
sion, respectively (see Table [l]). The gap between days 
35 and 50 is partly due to a major solar flare. Below we 
refer to these two data sets as early and late periods. 

Unlike the early period, when emergence of the 56 Ni 
lines strongly depends on the distribution of the radioac¬ 
tive Ni through the ejecta, for the late period the emis¬ 
sion in 56 Co lines is a generic prediction of all plausible 
models. Two energy bands optimal for detection of the 


SN signal in gamma-rays are clear from Fig PEI These 

t s^Co lir 


two bands, containing the most prominent 


lines, 


were used to generate images. The images w ere extracte d 


from SPI data from the late period as in (Churazov et 
aT7||20 141)). Namely, we vary the assumed position of the 


source and repeat the flux fitting procedure (see (3.1l for 
each position. The resulting images of the signal-to-noise 
ratio in the 835-870 and 1220-1270 keV energy bands are 
shown in Figj6] In both energy bands the highest peaks 
(4.7 and 4.3 a respectively) coincide well (within 0.3° ) 
with the SN2014J position, marked by a cross. 

The ISGRI spectra extracted at the known position 
of SN2014J for the early and late period are shown in 
Fig. Low energy (less than ~70 keV) part of the ex¬ 
tracted spectrum is likely contaminated by other sources 
in M82. 


4. MODELS 

4.1. A set of representative models 

For comparison with the INTEGRAL data we used 
a set of representative ID models (Table [2h, based on 
calculations of explosive nucleosynthesis models. To the 
first approximation, these models are characterized by 
the amount of radioactive nickel, total mass of the ejecta 
and the expansion velocity. Although current state-of- 
the-art si mulations of type la explosions can be done in 


3D (e.g., Seitenzahl et al. 

2013 

Fink et al. 

2014 

Moll 

et al. 

2014p, using these models would introc 

luce an ad- 


chtional viewing angle dependence. In order to avoid 
this extra degree of freedom and given that the overall 



















































5 


TABLE 1 
Data sets 


Set 

Dates 

Days since explosion 

Exposure", Msec 

early 

2014-01-31 : 

2014-02-20 

16 : 

35 

1.0 

late 

2014-03-05 : 

2014-06-25 

50 : 

162 

4.3 


“ Corrected for the periods of high background and the dead-time of SPI 



Fig. 2.— Estimated contribution of the off-diagonal terms in the 
SPI spectral response to the SN spectrum. The blue line shows the 
predicted spectrum of the W7 model for the late period, convolved 
with a simplified (nearly diagonal) SPI response. In this approxi¬ 
mation the instrumental broadening is parametrized as an energy 
dependent Gaussian with the width according to eq[2] The red 
line shows the same spectrum convolved with the response which 
includes estimated off-diagonal terms, caused by Compton scatter¬ 
ing of incident photons in the detector and surrounding structures. 
The off-diagonal component alone is shown with the dashed black 
line. The off-diagonal terms create a long low-energy tails associ¬ 
ated with gamma-ray lines. The impact on the brightest lines is 
small, while the continuum is strongly affected, especially at low 
energies. The model W7 is averaged over the period 50-162 days 
after the explosion. 


TABLE 2 

Set of models used in the paper 


Model 

MjVi, Mq 

Mf.ot, Mq 

E k , 10 51 erg 

DDTlpl 

0.54 

1.36 

1.29 

DDTlp4halo 

0.62 

1.55 

1.3 

DDTe 

0.51 

1.37 

1.09 

DETO 

1.16 

1.38 

1.44 

HED6 

0.26 

0.77 

0.72 

W7 

0.59 

1.38 

1.24 

ddtlp4 

0.66 

1.36 

1.35 

3Dbbal 

0.66+0.04“ 

1.36 

1.35 

DD4 

0.61 

1.39 

1.24 


“ additional “plume” of 56 Ni. 


significance of the SN2014J detection in gamma-rays by 
INTEGRAL (see fj5] and (J5| corresponds to only ~ 10 
s.t.d., we decided to keep m this work only a set of ID 
models to confront with the data. 


T he set of models in cludes the deflagration model 
W7 (Nomoto et a l. 1984), pure detonation model DETO 
(IBade nes et al . 2003), the sub -Chandrasekhar model 
HED6 (|iioefficFs Khokhlov|l996), an d several variants of 
the d elayed d etonation models: D D4 (Woosley & Weaver| 
1991|), DDTe (IBa denes et al.|2003|), DDTlpl, DDTlp4halo, 
ddtlp4, 3Dbbal (Isern et al. 20151. The ddtlp4 model 
was built to match the mass of b6 Ni suggested by the 
early optical evolutio n of SN2 014 J as d etected with the 
OMC of INTEGRAL (|Isem et al.|([2014)); P. Hofflich, pri¬ 
vate communication). In it, the transition density from 
deflagration to detonation was fixed at 1.4 10 7 g cm -3 . 
Model DDTlp4halo is a variant of the later in which 
the white dwarf is surrounded by a 0.2 Mq envelope, 
as might result from a delayed merger explosion. The 
3Dbbal model is essentially the same as the ddtlp4 plus 
a plume of 0.04 Mq of radioactiv e 56 Ni receding from 
the observer (see Isern et ah 2015 for details). 

The emerging X-ray ancT gamma-ray radiation from 
the expanding SNIa is determined by the total amount 
of radioactive isotopes, their distribution over velocities, 
the mass and the chemical composition of the ejecta and 
expansion rate. The processes are essent ially th e sam e as 
in type II supernovae (see, e 


et al. 


1987 


for 

a prototypical example of type 11 supernova - fcfiN19^7A). 
However, the mass of the ejecta and expansion rate differ 


.g. Sunyaev 
re il superr 


gamma-rays (see, e.g. 

Clayton 

et al. 

1969 

Woosley et al. 

1981 

Ambwani & Sutherland 

1988 

. A comprehensive 


different repres entative models was presented in The & 
Burrows (2014). 


Here we use the results of similar calculations (see be¬ 
low), which account for line broadening, needed for sys¬ 
tematic comparison with the INTEGRAL data. 

A Monte-Carlo code follows the propagation of the 
7—photons through the ejecta and accounts for scatter¬ 
ing and photoabsorption of photons and annihilation of 
positrons. The predicted spectra were generated with a 
time step of one day, covering the entire observational 
period. These model spectra were then averaged over 
the periods of 16-35 and 50-162 days respectively, to pro¬ 
vide fair comparison with the INTEGRAL results for the 
early and late periods. In particular, the effect of vary¬ 
ing opacity in each model over the observational period is 
correctly captured by this procedure. The computations 
include full treatment of Compton scattering (coherent 
and incoherent ), photoabsorption and pair production 
(see Milne et al. 2004 for details). The positrons pro¬ 
duced byyd^^ecayoT 36 Co (19% of all decays) annihilate 
in place via positronium formation. Both two-photon an¬ 
nihilation into the 511 keV line and the orthopositronium 
continuum are included. 


4.2. Transparent ejecta model (TEM ) 
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0 50 100 150 


Time, days 


Fig. 3.— ISGRI light curve in the 100-200 keV band. The S/N ratio in this band is expected to be the highest for the plausible models. 
The curves show the expected flux evolution for a set of models (see Color coding is explained in the legend. 


As we discuss below (p]) 
broadly consistent with the 
Table |2] However 


the INTEGRAL data are 
subset of models listed in 


Diehl et al. (2014) reported an evi¬ 


dence of 56 Ni at the surface in the first observation s of 
SN2014J with INTEGRAL (see also Isern et al.|2015| for 
an alternative analysis of early SN2014J observations). 
Presence of radioactive material at the surface would be 
an important result, since traditional models, listed in 
Table [2] do not predict it. One can attempt to patch 
our ID models with an additional component describing 
an extra radioactive material at the surface. Assuming 
that the material at the surface is transparent to gamma- 
rays, the fluxes of individual lines associated with Ni and 
Co decay, their energies and widths can be tied together. 
The transparency assumption is justified by the large ve¬ 
locities and small initial densities expected for matter at 
the surface of supernovae ejecta. In any case, it provides 
a lower limit to the mass of radioactive material, as opac¬ 
ity would demand a larger gamma-ray production rate in 
order to explain a given gamma-ray flux. This approach 
allows to describe many lines, associated with a trans¬ 
parent clump with only 3 parameters. Below we refere 
to this component as a Transparent Ejecta model (TEM), 
and use it in combination with the best-perf ormi ng W7 
model from out default set ID models (see |4.1[ ), i.e., 
the data are compared with the predictions of W7+TEM 
model. While this model by itself is not selfconsistent, it 


can be used to answer the following questions: 

• Once the predicted signal for the W7 model is re¬ 
moved from the observed spectra, do residuals re¬ 
semble a signal expected from a transparent clump 
of radioactive material? 

• Given the statistics accumulated by INTEGRAL , 
how much radioactive material in a transparent 
clump can be “hidden” in the data on top of a 
given ID model? 

In this section we describe the TEM model and then apply 


it to the data in 15.3 


The TEM model assumes that all line energies are shifted 
proportionally to their energies (i.e., the same velocity 
structure for all li nes), w hile their flux ratios follow the 
predicted ratios (Nadyozhin 1994) based on the decay 
chains of ^Ni-^Co-^Fe and OY Ni-> 57 Co^ 57 Fe. The 
list of the lines and their fluxes normalized to 1 M 0 of 
56 Ni are given in Table [3j For a given time period the 
model has 3 parameters: the initial 56 Ni mass (Mm), 
energy/redshift of the 847 keV line (E$ 47) and the broad¬ 
ening of the 847 keV line (0847)- The width of each line 
(Gaussian a) is defined as 


&line = (7847 x 


f Rline \ 
\ F/847 / 


(3) 
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Fig. 4.— The same as in Fig|3]for SPI data in two narrow bands near the brightest 56 Co lines. 


Ortho-positronium continuum and pair production by 
gamma-ray photons are neglected, while the 511 keV 
line is added assuming that 19% of 56 Co decays produce 
positrons, of which 25% form para-positronium yielding 
two 511 keV photons. 

5. RESULTS 

5.1. Combined ISGRI+SPI spectrum 
The SPI images (Fig|6j) for late period unambiguously 
show the characteristic signatures of 56 Co decay from 
SN2014J. A more quantitative statement on the amount 
of 56 Ni synthesized during explosion and on the proper¬ 
ties of the ejecta can be obtained from the comparison 
of the data with the predictions of the models. Since 


the late period is less affected by the transparency of 
the ejecta we start our analysis with the total spectrum 
obtained by INTEGRAL over this period. 

5.1.1. Late data 

The results of fitting of the combined ISGRI+SPI spec¬ 
trum (Figj8j) for the late period are given in Table [4] A 
full set of models from Tabj2]is used. The two groups of 
columns in Table [3] differ by the energy range in the SPI 
data used for comparison with the model. In the first 
group the data of ISGRI (70-600 keV) and SPI (400- 
1350 keV) are used. The data below 70 keV are likely 
contaminated by other sources in M82. The SPI data be¬ 
low 400 keV are omitted since during the late period the 
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Fig. 5. — early and late periods of INTEGRAL observations 
used for spectra extraction, shown as thick horizontal bars. Three 
curves show the evolution of the 56 Ni, 5 ®Co and 56 Fe masses, re¬ 
spectively, normalized to the initial 56 Ni mass. Note that opac¬ 
ity effects tend to suppress the emergence of gamma-rays at early 
phases of the supernova evolution, unless radioactive isotopes are 
present in the outer layers of the ejecta, or the explosion is strongly 
asymmetric. The dashed red line shows the 56 Co mass scaled down 
by the ratio of Co and Ni decay times TCo/rNi, which allows one 
to compare the expected relative strength of Ni (blue curve) and 
Co (dashed red curve) gamma-ray lines as a function of time. 


TABLE 3 

Line fluxes averaged over days 50-162 for a 

TRANSPARENT EJECTA MODEL (TEM ) FOR THE 
INITIAL 1 Mg OF 56 Ni 


f^Hne 1 keV 

Fline / ^847 

Flux“ 

Isotope 

846.78 

1.00 


6.57 

10" 4 

56 Co 

158.38 

7.98 

10“ 3 

5.25 

10“ 6 

56 Ni 

1561.80 

1.12 

10“ 3 

7.34 

10- 7 

56 Ni 

749.95 

3.99 

10 -3 

2.62 

10~ 6 

56 Ni 

269.50 

2.87 

10" 3 

1.89 

10~ 6 

56 Ni 

480.44 

2.87 

10 -3 

1.89 

10“ 6 

56 Ni 

811.85 

6.86 

10 -3 

4.51 

10“ 6 

56 Ni 

511.00 

9.50 

10 -2 

6.24 

10- 5 

56 Co 

1037.83 

1.40 

10" 1 

9.20 

10- 5 

56 Co 

1238.28 

6.80 

10- 1 

4.47 

10- 4 

56 Co 

*1771.49 

1.60 

10- 1 

1.05 

10- 4 

56 Co 

*2034.92 

7.90 

10“ 2 

5.19 

10- 5 

56 Co 

*2598.58 

1.69 

10- 1 

1.11 

10- 4 

56 Co 

*3253.60 

7.40 

10 -2 

4.86 

10- 5 

56 Co 

*14.41 

1.19 

10 -3 

7.80 

10- 7 

57 Co 

122.06 

1.03 

10 -2 

6.79 

10“ 6 

57 Co 

136.47 

1.19 

10" 3 

7.80 

10- 7 

57 Co 


Note. — “ - Flux is in units of phot s 1 cm 2 
* - Line is outside the energy range used for fitting 


data at these energies are expected to be dominated by 
the off-diagonal response of SPI. I.e. the observed SPI 
spectrum below 500 keV includes significant contribu¬ 
tion of the gamma-ray photons at higher energies, which 
are down-scattered inside the body of the telescope (see 
Fig(2|. The Null model (no source) gives y 2 = 1945.38 
for 1906 spectral bins. The improvement of the y 2 rela¬ 
tive to the Null is calculated by fixing the normalization 
at the predicted value for D = 3.5 Mpc (column 2) and 
by letting it free (columns 3 and 4). The typical value of 
the Ay 2 ~ 65 suggests ~ 8 a detection. 

One can draw two conclusions from this exercise. First 
of all a set of canonical ID deflagration (W7) or delayed 
detonation models (e.g., DD4) fit the data well without 
any adjustments to the normalization. The pure det¬ 
onation model DETD and a sub-Chandrasekhar model 
HED6 give poor fit and overproduce/underproduce the 
observed flux, respectively. Secondly, once the normal¬ 
ization is allowed to vary, all models give almost identi¬ 
cal gain in the y 2 , suggesting that relative strength of all 
prominent features is comparable in all models. Given 
the uncertainty in the distance to SN2014J (or M82) a 
deviation of the normalization at the level of ~20% can 
not be excluded. But DETO and HED6 models require by 
far larger changes in the normalization. 

While in the above analysis the SPI data with E < 400 
keV have been omitted to concentrate on the data less 
affected by the off-diagonal response, the right part of 
Tab|4] extends the analysis down to 70 keV for both in¬ 
struments. The basic conclusions remain the same, al¬ 
though, as expected, the significance of the detection in¬ 
creases to > 9 a. 


5.1.2. Early data 

We now proceed with the same analysis of the early 
data. Table [5] contains the gain in the y 2 for the same 
set of models. 

The DETO is clearly inconsistent with the data - inclu¬ 
sion of the model increases the y 2 relative to the Null 
model (no source). The HED6 model, which gave a poor 
fit to the late data, yields the y 2 comparable to other 
models. This is because smaller amount of 56 Ni is com¬ 
pensated by larger transparency of the lower-mass ejecta, 
which is important for the early data. 

The 3Dbbal model gives poor gain in y 2 if the nor¬ 
malization is fixed and the SPI data below 400 keV are 
excluded. If the normalization is free, and, especially, 
if the SPI data below 400 keV are included, this model 
performs marginally better than other models. However, 
it performs significantly better than other models when 
the SPI data below 400 keV are included. This is not 
surprising, since 3Dbbal model has been desi gned to fit 
the SPI data during this period (see|Isern et ah 2015, for 
details). The different “ranking” of the 3Dbbal model 
seen in Table [5] when SPI data below 400 keV are in¬ 
cluded or excluded, suggests a tension in the comparison 
of the fixed-normalization 3Dbbal model with the SPI 
and ISGRI data and also with the SPI data below and 
above 400 keV (see below). 


5.1.3. Early and Late data together 

Finally, in Table [6] we compare jointly the early and 
late data of ISGRI and SPI with the models, calculated 
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Fig. 6. — SPI images (S/N ratio) during late period in two narrow bands around most prominent 56 Co lines. Contours are at 2, 2.5 ... 5 
<j. Cross shows the position of SN2014J. The brightest peaks in each image coincide well with the position of SN2014J. Due to the dither 
patterrp] used during observations of SN2014J the central part of the image is much better covered than the outer regions. It is therefore 
not surprising that the level of noise is increasing away from the nominal target. 
a ht t p: / / www. cosmos. esa. int / web / integral 
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Fig. 7.— ISGRI spectrum measured at the position of SN2014J during early (red) and late (blue) periods. The energies of the second 
set of points are multiplied by a factor 1.02 for the sake of clarity. Dashed histograms show the predicted spectra of the W7 model for the 
same periods. The agreement with the predictio ns is reasonable excep t for the energies lower than 70 keV, where the spectrum is likely 
contaminated by other sources in M82 (see, e.g, |Sazonov et al.||2014| >. Dark green line shows crude approximation of the M82 spectrum 
measured before 2014. 
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Fig. 8 .— Combined ISGRI/SPI spectrum for the late period. The model (W7, see Tab[2|) has been convolved with the SPI off-diagonal 
response. The SPI data below 450 keV are omitted since during late period the data at these energies are expected to be dominated by 
the off-diagonal response of SPI. 


TABLE 4 

Ax 2 FOR BASIC MODELS FOR FIXED AND FREE NORMALIZATION RELATIVE TO THE NULL 
MODEL OF NO SOURCE FOR THE late PERIOD. 


Dataset: 

ISGRI(70-600 keV)&SPI(400-1350 keV) 

ISGRI(70-600 keV)&SPI(70-1350 keV) 

Model 

N = 1, Ax 2 

Nf 

1y jree 

Ax 2 

N = 1, Ax 2 

Nt 

1 ’free 

Ax 2 

DDTlpl 

66.4 

1.03± 0.13 

66.5 

87.3 

1.09± 0.12 

87.9 

DDTlp4halo 

65.9 

0.89± 0.11 

66.9 

88.1 

0.93=t 0.10 

88.5 

DDTe 

62.1 

1.09± 0.14 

62.5 

82.3 

1.15± 0.13 

83.7 

DETQ 

10.1 

0.52± 0.06 

66.4 

30.2 

0.55=t 0.06 

87.7 

HED6 

47.8 

1.86± 0.24 

60.7 

60.1 

2 .01± 0.22 

80.5 

W7 

65.0 

0.94± 0.12 

65.3 

86.9 

1 .01± 0.11 

86.9 

ddtlp4 

64.9 

0.85± 0.10 

66.9 

87.4 

0.90=t 0.10 

88.4 

3Dbbal 

63.2 

0.83± 0.10 

66.1 

85.7 

0.88± 0.09 

87.5 

DD4 

64.7 

0.89± 0.11 

65.7 

87.0 

0.95± 0.10 

87.3 

No source, x 2 (d.o.f.) 


1945.4 (1906) 



2696.9 (2566) 



Note. — AT is the normalization of the model with N = 1 corresponding to the explosion 
at the distance of 3.5 Mpc. Ax 2 characterizes an improvement of x 2 f° r a given model 
relative to the Null model. Larger positive values indicate that the model is describing 
the data significantly better than other models (see Appendix). The data below 70 keV are 
likely contaminated by other sources in M82. SPI data below 400 keV are omitted in the first 
dataset (left half of the Table) since the data at these energies are expected to be dominated 
by the off-diagonal response of SPI (see j |3.1[ ). 
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Energy [keV] 


Fig. 9. — Combined ISGRI/SPI spectrum for the early period. The model (W7, see Tab[2]l has been convolved with the SPI off-diagonal 
response. 


TABLE 5 

The same as in Table [4] for the early period. 


Dataset: ISGRI(70-600 keV)&SPI(400-1350 keV) ISGRI(70-600 keV)&SPI(70-1350 keV) 


Model 

N = 1, Ay 2 

1v free 

Ay 2 

N = 1, Ax 2 

1 ’free 

Ax 2 

DDTlpl 

14.9 

0.84± 0.21 

15.4 

33.2 

l.lli 0.19 

33.5 

DDTlp4halo 

14.6 

1.00± 0.26 

14.6 

29.8 

1.34± 0.24 

31.8 

DDTe 

14.3 

1.30± 0.33 

15.1 

26.9 

1.72± 0.30 

32.6 

DETD 

-83.9 

0.28± 0.07 

14.8 

-64.8 

0.37± 0.06 

35.2 

HED6 

15.7 

1.05± 0.26 

15.8 

32.7 

1.39± 0.23 

35.5 

W7 

15.9 

0.87± 0.22 

16.2 

34.8 

1.14± 0.19 

35.3 

ddtlp4 

11.3 

0.65± 0.17 

15.7 

33.3 

0.86± 0.15 

34.2 

3Dbball 

6.7 

0.56± 0.13 

17.6 

37.0 

0.76± 0.12 

41.4 

DD4 

14.1 

0.77± 0.19 

15.5 

33.6 

1.01± 0.17 

33.6 

No source, x 2 (d.o.f.) 


1856.7 (1906) 



2615.9 (2566) 



for corresponding periods. The two columns in Table [6] 
differ by the energy range in the SPI data used for com¬ 
parison with the model. In each case the normalization 
was fixed at the value set by the adopted distance of 3.5 
Mpc. In each column we mark with bold face the mod¬ 
els which have Ay 2 different from the model with the 
largest Ay 2 by less than 4 (see Appendix for the clarifi¬ 
cation on the interpretation of this criterion in Bayesian 
and frequentist approaches). Once again, ID deflagra¬ 
tion model W7 and ” standard” delayed detonation model 
perform well. The 3Dbbal, which was designed to ac¬ 
count for tentative feature in the early SPI data at low 
energies, not surprisingly perforins well if the SPI data 
below 400 keV are included. However, if only the data 


above 400 keV are used for SPI, this model yields signif¬ 
icantly lower Ay 2 than the W7 or DDTlpl models. 

5.2. Comparison of gamma-ray light curves with models 

While the spectra for the early and late periods al¬ 
ready provide an overall test of the basic models, addi¬ 
tional information can be obtained by analyzing the time 
variations of the fluxes in broad energy bands (see Figj3] 
and [4]). The total number of time bins is 34. Each bin 
corresponds to one revolution (i.e. ~3 days). The first 
raw in Table [7] provides the values of the y 2 (for Null 
model of no source) in three energy bands: 100-200 keV 
(ISGRI), 835-870 keV (SPI) and 1220-1272 keV (SPI). 
The normalization of the model lightcurves is fixed to 1. 
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TABLE 6 

Ax 2 FOR THE JOINT DATA SET OF THE early AND late SPECTRA FOR A BASIC SET OF 
MODELS WITH FIXED NORMALIZATION. THE VALUE OF Ay 2 SHOWS THE IMPROVEMENT OF THE 
X 2 RELATIVE TO THE NULL MODEL OF NO SOURCE. 


Model 

ISGRI & SPI(400-1350 keV) 

Ax 2 

ISGRI & SPI(70-1350 keV) 
Ax 2 

DDTlpl 

81.3 

120.5 

DDTlp4halo 

80.5 

117.8 

DDTe 

76.4 

109.2 

DETO 

-73.8 

-34.7 

HED6 

63.5 

92.8 

W7 

80.9 

121.7 

ddtlp4 

76.2 

120.7 

3Dbball 

69.9 

122.7 

DD4 

78.8 

120.7 


Note. — Bold-faced are the models which have Ax 2 different from the model with the 
largest Ax 2 by less than 4, the criterion used to group models into “more plausible” and 
“less plausible” respectively (see Appendix). 



Fig. 10.— Spectra predicted by the TEM model for early (blue) 
and late (red) data sets, convolved with the SPI response. The 
broadening of the reference 847 keV line is set to 20 keV (Gaussian 
sigma). The initial 56 Ni mass is 1 Mq. 

For 34 bins the value y 2 for a correct model is expected 
to be in the interval ~26-42 in 68% of cases. Clearly, the 
Null model does not fit the data well. 

Other raws show the improvement of the y 2 relative to 
the Null model. I.e., Ay 2 = X NuU - xL del- From Ta " 
ble [7] it is clear that DETO model strongly overpredicts the 
flux in all bands and can be excluded (y 2 becomes worse 
when this model is used). Other models leads to signifi¬ 
cant improvement with respect to the Null model, except 
for the 3Dbbal model in the 100-200 keV band where it 
exceeds the observed flux in the early observation, while 
in the SPI bands all these models are comparable. 

The last column in Table [7] provides the y 2 for three 
bands joints. This is basically the sum of the values of y 2 
for individual bands. Bold-face d are the best performing 
models: W7 and DDT1P1. As in §5. 1| these are the models 
which have Ay 2 different from the model with the largest 


Ay 2 by less than 4 (see Appendix). 

One can also compare the lightcurves with the hypoth¬ 
esis of a constant flux. The mean level of flux was esti¬ 
mated for each band and the value of the y 2 was cal¬ 
culated. The values of Ay 2 relative to “No source” are 
given in the last row of Table [7] One can see that this 
simple model is almost as good as other best-performing 
models in individual bands (even taking into account 
that this model has a free parameter - mean flux). This 
is of course the result of low statistical significance of the 
SN2014J detection that makes it difficult to constrain 
time variations of a faint signal. For the combined values 
for all three band the effective number of free parameter 
is 3 (mean fluxes in each band) and one can conclude 
that, e.g. W7 model performs marginally better than the 
constant flux model. 


5.3. Search for the velocity substructure in the late data 

The above analysis suggests that the INTEGRAL data 
broadly agree with a subset of simple ID models (e.g., 
W7 or DD4). Since the true structure of SN2014J is surely 
more complicated than predicted by ID models, it is in¬ 
teresting to verify if adding an extra component to the 
model (on top of the best-performing W7 model) signifi¬ 
cantly improves the fit. In this section we use TEM model 
as such extra component. This choice is partly driven by 



lines associated with the Ni—>Co—>-Fe decay in the TEM 
model are tied to the energy (redshift) and the width of 
the reference 847 keV line. The flux ratios are also tied 
together using a model of an optically thin clump, tak¬ 
ing into account time evolution of the Ni and Co masses. 
Examples of spectra predicted by TEM model (for 1 Mq 
of 56 Ni) are shown in Fig. 


10 


Thus, we consider a composite model, consisting from 
the W7 model (with the normalization fixed to 1) and the 
TEM model. This two-component (W7+TEM ) model ef¬ 
fectively searches for a transparent clump of radioa ctiv e 
material on top of the base-line W7 model (see FigpTj). 
The horizontal axis shows the energy of the reference 
847 keV line in the observer frame and different colors 
correspond to different 847 keV line broadening param- 
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TABLE 7 

A,\ 2 FOR LIGHT-CURVES IN THREE ENERGY BANDS FOR DIFFERENT MODELS. THE VALUE OF A\' 2 SHOWS THE 
IMPROVEMENT OF THE X 2 RELATIVE TO THE NULL MODEL OF NO SOURCE. THE VALUE OF THE \ 2 FOR THE 

Null model is given in the first raw. 


Model 

100-200 keV (ISGRI) 

835-870 keV (SPI) 

1220-1272 keV (SPI) 

Three bands jointly 

No source, y 2 

51.0 

49.6 

51.1 

151.7 

DDT1P1 

16.9 

18.6 

19.3 

54.8 

DDTlP4halo 

9.2 

17.4 

19.5 

46.1 

DDTe 

16.7 

17.1 

17.6 

51.4 

DET0 

-105.0 

-6.3 

13.5 

-97.8 

HED6 

20.6 

16.0 

14.1 

50.7 

W7 

18.8 

18.4 

20.0 

57.2 

DDT1P4 

9.1 

17.5 

20.3 

46.9 

3Dbball 

-4.8 

17.4 

20.4 

33.0 

DD4 

12.7 

17.9 

20.0 

50.6 

CONST 

18.7 

16.1 

18.5 

53.3 


Note. — Total number of time bins is 34. For the joint x 2 the effective number of bins is three times 
larger - 102. 


eterized through a Gaussian a - see legend. For a given 
redshift/energy and width of the reference 847 keV line 
the model has the normalization (initial 56 Ni mass) as 
the only free parameter. The best-fitting 56 Ni mass is 
sho wn i n the top panel of FigJTT] The bottom panel 
(Fig |ll| ) shows the improvement in the y 2 (relative to 
the W7 model alone) due to the TEM model. 

As is clear from Fig. [ll] this model does no provide 
compelling evidence for a transparent clump on top of 
the W7 in the late data. Formally, there is a Ay 2 ~ 9.5 
peak at ~858.5 keV, which corresponds to a narrow (~1 
keV broad, red curve) component with a negative mass 
of —0.05 Mq, which can be interpreted as a marginal 
evidence for a dip in the velocity substructure, given 
that this improvement of the Ay 2 came at the cost of 
adding three more parameter^] to the model. One can 
estimate the constraints on the line flux (mass of a trans¬ 
parent clump) that such analysis can provide, by fix¬ 
ing the centroid energy and the width of the reference 
847 keV line and calculating the expected statistical un¬ 
certainty. Since the normalization of the TEM model is 
the only free parameter in this particular experiment, 
the estimation of the uncertainty is straightforward (see 
Fig [l2]) • Three curves shown in Fig[l2] show lcr uncer¬ 
tainty on the initial 56 Ni mass for the early set (dashed- 
blue: SPI data in the 70-1350 keV band; long-dashed- 
green: 400-1350 keV) and late set (solid-red: 400-1350 
keV), respectively. Conservative upper limit based on 
the assumption of pure statistical errors would be 3 times 
these values. Letting the broadening and the redshift to 
be free parameters (look-elsewhere effect) would increase 
this limit even further. 

These experiments show that the late data are consis¬ 
tent with a presence of a velocity substructure (parame¬ 
terized via our TEM model) on top of the ID W7 model at 
the level ~ 0.05 Mq, provided that the lines are slightly 
broadenened. 

We now do a similar experiment with the early data, 
using a TEM +W7 model for SPI data in the 400-1350 keV 
band (Fig|l3| left panel) and in the 70-1350 keV band 
(Fig[l3j right panel), respectively. 

3 We note, that the width and especially energy of the reference 
line are very nonlinear parameter that could lead to large changes 
in the y 2 . 



750 800 850 900 950 


E, keV 


Fig. 11. — Fitting the SPI data in the 400-1350 keV band with 
a composite W7 + TEM model. The normalization of the W7 model is 
fixed to 1. In the TEM model all lines are tied to the energy (redshift) 
and the width of the reference 847 keV line. The flux ratios are 
tied using a model of optically thin clump, taking into account 
time evolution of the Ni and Co masses. This setup is optimized 
for a search of a transparent clump on top of the W7 model). For a 
given energy and width the model has only normalization (initial 
56 Ni mass in the clump) as a free parameter. The bottom panel 
shows the improvement in y 2 and the top panel shows the best¬ 
fitting 56 Ni clump mass. Different colors correspond to a different 
847 keV line broadening parameterized through a Gaussian a - see 
legend. No compelling evidence for a clump is seen in the data. The 
sensitivity of the data to the mass of the clump strongly depends 
on the broadening of the lines (see Fig |12| ). 

The left panel does not show any significance evidence 
for a clump on top of the W7 model. The structure in the 
right panel is more complicated. The data used in this 
panel now include the 56 Ni line at 158 keV. We note, that 
if the 158 keV line is able to escape, then it is certainly 
true for higher energy lines of Ni and Co. Therefore the 
analysis should be done for the whole band to achieve 
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Fig. 12. — Uncertainty in the initial 56 Ni mass as a function of 
line broadening for early set (dashed-blue: SPI data in the 70- 
1350 keV band; long-dashed-green: 400-1350 keV) and late set 
(solid-red: 400-1350 keV), respectively, assuming transparency to 
gamma-rays generated close to the surface. A conservative up¬ 
per limit on the initial mass of ’’extra” radioactive 56 Ni, is three 
times this value at a given line width. For the line broadening of 
10 4 km s _1 (FWHM), the expected value of (T 847 is ~ 12 keV. This 
value can be regarded as a fiducial value for a simple SNIa model. 


the most significant results. First of all, our analysis 
does not show compelling evid ence for a n arrow and un¬ 
shifted component reported in Diehl et al.j (20141 - there 
is a weak (Ay 2 ~ 6, i.e. ~ 2.4 a detection if we ignore 
the freedom in the redsliift and broadening) peak at 847.5 
keV, corresponding to a narrow line (black curve) with 
a mass of ~ 0.027 M 0 of 56 Ni. There are several sep¬ 
arate peaks of similar magnitudes, covering the energy 
range of interest. However there is a more significant (al¬ 
beit also marginal) evidence for a redshifted and broad 
component with Mm ~ 0.0 8 M 0 , E ~ 826.5 keV and 
a ~ 8 keV (see Isern et al. 2015 for discussion). The 
gain in y 2 is ~l8 and for a fixed energy and broadening 
(putting under the rug possible systematic errors in the 
background modeling and uncertainties in the calibration 
of the off-diagonal response) this would be a 4.2 a detec¬ 
tion. However the freedom in the energy, width (look 
elsewhere effect) and the normalization deteriorates the 
significance. Should all these free parameters be linear 
(as is normalization), one would expect the change in the 
y 2 of ~ 3 due to pure statistical fluctuations. However, 
the energy and the width are nonlinea r an d t he g ain in 
y 2 might be significantly larger. In Fig. 11 and 13 we see 
multiple peaks with the change/gain in y 2 up to ~10. 
Assuming that the latter value can be used as a crude 
estimate of a possible gain in y 2 due to non-linearity of 
the TEM model, the significance of the detection of the 
excess drops below 3<r. 

Taking the best-fitting parameters at the face value, we 
can go back to the late data and compare the spectra (in 


the 400-1350 keV band) with the TEM +W7 model, freezing 
TEM model parameters at the best-fitting values obtained 
for the early data. This gives the y 2 = 1883.05, i.e. 
worse than the W7 model alone (y 2 = 1879.3). If we let 
the normalizations of both TEM and W7 models free (but 
freezing energy and broadening of the TEM model), then 
we improve slightly the y 2 to 1878.9, but the best-fitting 
mass becomes slightly negative, although consistent with 
zero —3 10 -3 ± 5 10 -2 M 0 , while the best-fitting nor¬ 
malization of W7 model becomes 0.92 (c.f. Tab. El where 
SPI data are used together with the ISGRI dataj! 

We concluded that there is a tension between “low” en¬ 
ergy SPI data in early observations and the rest of the 
INTEGRAL data (Tab. J5] and Tab. 0- However, this 
tension is not prohibitively large and could be attributed 
to statistical fluctuations in the data, if a conservative 
approach is adopted. A possible evidence of the red- 
shifted and broadened 158 keV line in the early data 
and p ossib le implications are further discussed in |Isern| 


et al. (2015). 


5.4. SPAR model 


Apart from the models discussed above , we als o used 
3PAR model, introduced in Churazov et al. (2014b). This 
is a spherically symmetric model of homologously ex¬ 
panding ejecta with exponential density profile p ex e v ^ Vc . 
The model is cliaaracterized by three parameters: initial 
mass of the 56 Ni M Ni , total mass of the ejecta M e j ecta , 
and characteristic expansion velocity V e in the exponen¬ 
tial density distribution. In this model a mass-weighted 


root-mean-squared velocity of the ejecta is y/l2V e . 

The main shortcoming of this model is the assump¬ 
tion that all elements, including radioactive Ni and Co, 
are uniformly mixed through the entire ejecta. This is 
an ad-hoc assumption, made in order to stay with only 
three-parameteres model, but it is not justified. It has 
the major impact for the early gamma-ray light curve, 
producing gamma-ray emission even at the very early 
phase (see Fig. [TJ|3| [4b . At later times (day 50 or later), 
the role of mixing isless significant. We therefore ap¬ 
plied this model to the late ISGRI and SPI spectra to 
get estimates of Mm, M e j ecta and V e , which are not lim¬ 
ited to values characteristic to the set of plausible models 
given in Table [2j The main purpose of using this model 
is to understana the level of constraints provided by the 
INTEGRAL data on the main characteristics of the su¬ 
pernova. Simplicity of the model allows us to calculate 
this model on a large grid of possible values of Mm, 
M e j ecta and V e - 

A Monte Carlo radiative transfer code is used to cal¬ 
culate the emergent spectrum, which includes full treat¬ 
ment of Compton scattering (coherent and incoherent) 
and photoabsorption. Pair production by 7-ray photons 
is neglected. The positrons produced by /3 + decay an¬ 
nihilate in place via positronium formation. Both two- 
photon annihilation into the 511 keV line and the ortho- 
positronium continuum are includ ed. 

The results are shown in Fig. [l4j The best-fitting 
values M m = 0.63 M 0 , M ejecta = 1.8 M 0 , V e = 
3 10 3 km s -1 are marked with a cross. The ler con¬ 
fidence contours (corresponding to Ay 2 = 1, i.e. for 
single parameter of interest) are shown with the thick 
solid line. Clearly, the Ni mass M^i and the character- 
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Fig. 13. — Same as in Fig |ll| but for the early SPI spectrum. LEFT: TEM +W7 model and the SPI data in the 400-1350 keV band. The 
normalization of the W7 model is fixed to 1. RIGHT: TEM +W7 model and the SPI data in the 70-1350 keV band. The normalization of 
the W7 model is fixed to 1. The low energy part of the SPI spectrum is included to make sure that the se Ni line at 158 keV is within 
the energy range probed. There is a marginal evidence of a redshifted (by ~8000 km/s) component with the width of ~ 8 keV (Gaussian 
sigma), corresponding to Mjvi ~ 0.08 Mq. See text for the discussion. 


istic expansion velocity V e are better constrained than 
the ejecta mass. This is not surprizing, given that the 
data averaged over the period 50-162 days after explosion 
are used, when the ejecta are relatively transparent for 
gamma-rays. As a results the flux in the lines depends 
primarily on the Ni mass, line broadening is set by the ex¬ 
pansion velocity, while ejecta mass influence mostly the 
amplitude of the scattered component, which declines 
with time relative to the ortho-positronium continuum 
when the optical depth declines. If we fix the poorly 
constrained ejecta mass to M e j ecta = 1.4 Mq, then the 
derived Ni mass is constrained to the range 0.54-0.67 Mq. 

For the set of models listed in the Table [2] we can csti- 

/ Ex 

mate the effective V e using the relation V e = J --—, 

V 6M e j ec ta 

valid for pure exponential model. The values V e vary 
between ~ 2580 km s -1 for DDTe to ~ 2960 km s -1 
for DET0 models and is equal to 2740 km s -1 and 
2820 km s -1 for W7 and DDTlpl respectively. Not sur¬ 
prisingly all “successful” models (e.g. W7 and DDTlpl) 
have their char acte ristic parameters well inside contours 
plotted in Fig. [l4j while DET0 and HED6 are far outside 
the contours, primarily because of Ni mass. 

5.5. Summary of model fitting 

The comparison of the INTEGRAL data with the sub¬ 
set of models (see the sections above) allows one to 
crudely rank the models according to their success in 
different tests. For each test (data set) we can choose 
the “best” model, which provides the largest improve¬ 
ment Ay 2 compared to Null model (or having the small¬ 
est y 2 for the lightcurves). We can then adopt an ad 
hoc definition that other models that have y 2 different 
from the best model by 4 (i.e. ~ 2 a confidence) are 


classified as “good”. Similar approach can be applied to 
the lightcurves in each band (Tab. [7]), by adding 4 to 
the minimal value of the y 2 among models. Applying 
this test to Tables 0-0 we conclude that W7 and DDTlpl 
pass all these tests, closely followed by DD4, ddtlp4, and 
then by DDTlp4halo and 3Dbbal. DET0 and HED6 fail 
most of the tests. Of course, given the uncertainties in 
the distance, background modeling and calibration is¬ 
sues, we can not reject models other than DET0 and HED6. 
E.g., if we let the normalization to be a free parameter 
(equivalent of a statement that the distance is highly 
uncertain) then most of the models become barely dis¬ 
tinguishable. We rather state, that a whole class of near- 
Chandrasekhar models provides a reasonable description 
of the data, with the W7 and DDTlpl being the most suc¬ 
cessful, closely followed by a broader group of delayed- 
detonation models. 

6 . CONSISTENCY WITH OPTICAL DATA 

We now make several basic consistency checks of 
gamma-ray and optical data, using optical observations 
taken quasi-simultaneously with INTEGRAL observa¬ 
tions. 


6.1. Optical and gamma-ray luminosities 
We use BVRIJHK photometry reported by 


Foley et| 


unnnos- 


al. |2014) to estimate the bolometric (UVOIR) 
lty of SN 2014J on days 73 and 96 after the explosion. 
Since the data do not contain the (7-band photometry, 
we include the U magnitude recovered on the bases of the 
U — B color of the dereddened normal SN la, SN 2003hv 
(Leloudas et al. 2009). The SN 2014J fluxes were cor¬ 
rected for the extincti on using slightly differen t extinc- 
tion laws report ed by Amanullah et al. (|2014) and Fo¬ 
ley et al. (2014). The average of both fluxes for each 
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Fig. 14.— Confidence contours for 3PAR model, corresponding to A\ 2 = 1 with respect to the best-fitting value. The cross show the 
best-fitting parameters of the 3PAR model: My,; ~ 0.63 Mg, v e ~ 3000 km s -1 , M e j ecta ~ 1.8 Mg. The late ISGRI and SPI spectra 
are used for this analysis. Confidence intervals plotted in this figure correspond to 1 <r for a single parameter of interest. The largest 
uncertainty is in the mass of the ejecta, while the Ni mass is the best determined quantity. 


epoch were used then to produce the integrated flux. To 
this end we approximated the spectral energy distribu¬ 
tion by the combination of two functions each of which 
is a smooth broken power law. The SED integration in 
the range of 0.1 < A < 10 /tm with the distance of 3.5 
Mpc, results in the luminosity estimates of (11 ±1) x 10 41 
erg s _1 on day 73, and (6.5 ± 0.6) x 10 41 erg s _1 on day 
96. These values agree well with the estimated amount 
of deposited energy in the best-fitting 3PAR model: ~ 
1.0 x 10 42 erg s _1 and ~ 5.3 x 10 41 erg s _1 for day 73 
and 96 respectively. According to this model the fraction 
of thermalized energy is ~34% and ~20% for these dates 
respectively. 


6.2. Asymmetry in late optical spectra? 

The issue of asymmetry of SN 2014J ejecta is of vital 
importance because the strong deviation of the 56 Ni dis¬ 
tribution from the spherical symmetry would affect the 
interpretation of the ganrma-ray data. Generally, the 
asymmetry of the 56 Ni dist ribution is expected in the 
binary WD merger scenario (Pakmor et al. :2012|. More¬ 
over, a single degenerate scenario also does not rule out 
the ejecta asymmetry caused by the noncentral early de¬ 
flagration (Malone et al. 2014]) . In fact, signatures of 
asymmetry have been already detected in several SNe la 


at the nebular stage (t > 100 d). The asymmetry is man¬ 
ifested in the emissio n line shift or/and the double peak 
emission line profiles dM otohara et al.|2006 Maeda et al.| 
2010 |Dong et aL||2014b 


To probe a possible asymmetry of SN 2014J ejecta we 
rely on the nebular optical spectrum taken on day 119 
after the B maximum, i.e., 136 d after the explosion (Bik- 
maev et al. J2015) at the 1.5-m Russian-Turkish telescope 


(RTT-150) of the TUBITAK National Observatory (An¬ 
talya, Turkey). The SN 2014J spectrum corrected for 
the interstellar reddening in M82 of E(B — V) = 1 (c.f. 
Foley et al. 2014) is shown in Fig. |l5| together with that 


of SIN 201 lie obtained at the same instrument on day 141 
after the maximum. The spectra of both supernovae look 
similar except for the blueshift of SN 201 lfe emissions by 
~ 10 3 km s _1 relative to SN 2014J. 

We focus on the [Co III] 5890 A emission that is not 
hampered markedly by the blending with other lines. It 
should be emphasised that on day 136 d after the explo¬ 
sion this line is dominated by 56 Co; the contribution of 
57 Co and stable Co isotopes is negligible. The Thom¬ 
son optical depth at this epoch is small (~ 0.2) and does 
not affect the line profile. The [Co III] emission is the 
superposition of five lines of the a 4 F - a 2 G multiplet. 
Each line we describe by the Gaussian with the ampli- 
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Fig. 15.— The spectra of SN 2014J (day 119 after the maximum) 
and SN2011fe (day 14 1 after the maximum) obtained with RTT- 
150 telescope |Bikmaev et al. (|2015|). Overall the spectra are very 
similar in terms ot the flux level, line shape and line ratios. The 
exception is the prominent blueshift of [FeUI], [Fell], and [CoIII] 
emissions of SN 201 lfe relative to SN 2014J. The strong interstellar 
Nal absorption in the SN 2014J spectrum arises in the M 82 galaxy. 


tude proportional to the collisional excitation rate times 
the radiative branching ratio. We adopt the heliocentric 
recession velocity of +104 ± 15 km s -1 that takes into 
account the recession velocity of +203 km s -1 for M 82 
(NASA Extragalactic Database NED) and the rotational 
velocity of M 82 at the SN 2014J position. The best fit 
(Fig. 16) is found for the full width at half maximum 
for each line FWHM = 10450 km s -1 and the line shift 
of v s = +130 ± 17 km s -1 . With the exception of this 
small shift, each [Co III] line is fairly symmetric at least 
in the radial velocity range of |u r | < 6100 km s -1 . The 
small line shift may be related to either intrinsically small 
asymmetry of 56 Ni distribution, or the special viewing 
angle, if the ejecta is actually non-spherical. To sum¬ 
marize, the SN 2014J optical spectrum does not show 
signatures of strong asymmetry. 


7. DISCUSSION AND CONCLUSIONS 

We have analyzed a complete set of INTEGRAL ob- 
servations of SN2014J. We confirm our previous results 
(Churazov et al.|2014b) that the data are broadly consis¬ 
tent with the predictions of a nearly-Chandrasekhar WD 
explosion, with (ID) deflagration or delayed detonation 
models providing equally good description (see Tables [6] 
- [7| . While pure deflagration models are disfavored be¬ 
cause of the expected large scale mixing and incomplete 


Fig. 16.— [Colli] 5900 A emission in the SN2014J spectrum 
on day 119 (thin line) along with the model (thick line) which 
includes five component of the a 4 F - a 2 G multiplet. The narrow 
absorption feature at the top of the profile is due to Na I interstellar 
absorption in M 82. At the bottom shown is the residual ’’model 
minus observation”, which demonstrates a good fit in the range of 
5770-6060 A. 

burning in 3D simulations, in the ID case they yield the 
same gamma-ray flux as the delayed detonation models. 
Pure detonation (or strongly sub-Chandrasekhar) models 
strongly overproduce (underproduce) observed gamma- 
ray flux and can be excluded. Allowing a freedom in 
the normalization of the model (equivalent to allowing 
the initial mass of 56 Ni to be a free parameter, while 
keeping other parameters unchanged) makes all models 
essentially indistinguishable at the level of statistics, ac¬ 
cumulated by INTEGRAL . 

We have searched for possible velocity substructure on 
top of the predictions from ID models, by adding a set 
of broadened Gaussian lines to the best-performing W7 
model. The energies and fluxes of the lines are tied to 
the predictions of the Ni and Co decay chains, appropri¬ 
ate for the optically thin clump of Ni. This analysis did 
not reveal strong evidence for a prominent velocity sub¬ 
structure in the gamma-ray data during the late phase of 
the SN evolution (after day 50). Given the statistics ac¬ 
cumulated by INTEGRAL , a clump with the 56 Ni mass 
~ 0.05 Mq producing slightly broadened lines (FigJT2|) 
could be consistent with the late gamma-ray data. Sim¬ 
ilar analysis of the early data has a best-fitting solu¬ 
tion with a redshifted and broadened component with 
Mjvj ~ 0.08 Mq, E ~ 826.5 keV and a ~ 8 keV. How¬ 
ever, the statistical significance of this extra component 
is marginal and the late observations do not provide 
further evidence for the presence of such component (see 
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also |Diehl et aLl|2014| |Isern et ahj|2015[ for independent 
analysis of early observations of SN2014J). 

From the optical light curves and spectra SN2014J ap¬ 
pears to be a “normal” SNIa with layered structure and 


no evidence for large-sca le mixing (e.g., Marion et al. 


2015 Ashall et al. 


2014), consistent with the delayed- 


dctonation models . The detection of stable Ni (Friesen 


et al. 2014 Telesco et al. 20151 in IR suggests high 


density of the burning material, characteristic for near- 
Chandrasekhar WD. 

Optical spectrum taken at the nebular stage (day 
~ 136 after the explosion) also do not show strong asym¬ 
metry in the Co and Fe lines. Unless the viewing angle is 
special, the distribution of these elements in the ejecta is 
symmetric. These data do not provide any direct support 
for collision/merger scenario. The late SN2014J spec¬ 
trum is very similar to that of SN2011fe, albeit with the 
pronounced blueshift of emission lines of the latter. 

Apart from the above mentioned feature in the early 
observation, which we consider as marginal, the rest of 
the INTEGRAL and optical data appear consistent with 
the predictions of “canonical” ID explosion models of a 
nearly-Chandrasekhar carbon-oxygen white dwarf. 
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This work was based on observations with INTE- 


APPENDIX 

COMPARISON OF SIMPLE MODELS VIA A\ 2 CRITERION 

In this appendix we clarify our approach of comparing different models using ungrouped spectral and/or timing 
data. 


Grouping the data and small number of counts per bin 

Consider N data points A, i = 1, N and JV>1 (e.g., a spectrum measured in N energy bins) that correspond to 
a model M t , Di = M ti , + m , where is the noise with Gaussian distribution with zero mean and known variance 
erf. The noise is uncorrelated, i.e. < UiUj >= 0, if i / j. For simplicity we drop below the index i in the expressions 
containing summation over i. 

Suppose that we want to compare two competing models Mi and M 2 with no free parameters (this corresponds 
to Neyman-Pearson lemma of two simple hypotheses). For our purpose it is useful to write explicitly the probability 
distribution of Ay 2 between competing models. 

Let us define (5i = Mi — M t and S 2 = M 2 — M t and calculate y 2 for both models. For the model Mi 


Xi 


E 


(Mi — D 


= E 
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and a similar expression for M2. Thus Ax 2 = Xi ~ X 2 
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where <5 i j2 = Mi — M 2 . The last term in the above expression obviously has a Gaussian distribution with zero mean 


and variance 4^^ 


S\ 2 <n 2 > 
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. Thus 



(A4) 


where X has a Normal distribution. It is clear that the above expression does not depend on data grouping (see 
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Churazov et al., 1996) as long as the grouping does not severely affect the shapes of Si, 62 or <5^2. Furthermore, Ay 2 
can have a Gaussian distribution by the central limit theorem, even when the noise in the original data is not Gaussian. 
For instance, if n corresponds to the Poisson noise in the data with small number of counts per bin, the Ay 2 will 
have Gaussian distribution provided that the total number of counts contributing to E is large (Churazov et al., 
1996). We therefore can use the original data with no grouping to calculate Ay 2 . 


Dividing the models into a “more plausible” and “less plausible” groups 
For a given observed Ay 2 bs the ratio of likelihoods for two models is A = = e~ Ax ° bs ^ 2 , which corresponds 

2) 

to the Akaike information criterion (AIC) or Bayesian information criterion (BIC) criteria. Alternatively we can also 
employ Neyman-Pearson lemma to differentiate between two models. If M\ is the true model, then Mi = M t . <5i = 0, 
62 = —Si } 2 and the distribution of Ay 2 = y 2 — y 2 is 


Ay 2 = 


Then the probability a of getting Ay 2 > y is 



a = Q{x), 


(A5) 

(A6) 


where Q(x) = 1/2 erfc(x/-\/2), x 


A + y 

2VA 


and A = £ (¥)’•■ 


Although the value of A is known, we can take a 


conservative approach and write that > ^fy for any A E [0,00]. This value is achieved at A = y. Thus one can 
conclude that the conservative estimate of the probability of getting better y 2 for the wrong model (i.e., Ay 2 bs > y) 
corresponds to more than v fy standard deviations. In the paper we use the value of Ay 2 bs < 4 to separate the models 
into a “more plausible” and “less plausible” groups. As we emphasized above, one can also interpret this value in the 
frame of BIC or AIC approaches. 


Pairwise model comparison and the goodness of fit criterion 

When comparing two models in terms of the y 2 we are effectively projecting the N-dimensional data D on a 
line connecting the models Mi and M2. We can formally introduce a linear parameter p , make a composite model 
Mi + p(M 2 — All) and calculate best-fitting value pbf that minimizes the y 2 


Pbf 


E ~(D - Mi) 6 U 2 /a 2 
E<5iV^ 2 


(A7) 


It is obvious that this quantity has a Gaussian distribution with the standard deviation 

* P = (ZXa/* 2 ) -172 - ( A8 ) 

We can now use the deviation of pbf from 0 or 1 as a goodness of fit criterion for models Mi and M 2 , respectively. A 
significant deviation (i.e. \p b f/cr p \ 1 or \{p b f — l)/a p \ 1) implies that one of the models (or both) is unlikely. While 
the goodness of fit can be calculated for each model alone, the power of the criterion depends on the data grouping 
(imagine, for example, a very weak signal distributed over a large number of bins). The goodness of fit calculation 
described above verifies only the projection of the data on the line set by models, but it “optimally” compares the 
difference between plausible models with the signal present in the data. 

For each model used in 93 we have calculated the maximum deviation p b f /<j p | with respect to all other models 
(except for the 3Dbbal model that was designed to fit the INTEGRAL data). Table [S] provides corresponding values, 
when ISGRI and SPI data for the early and late periods are considered jointly. Based on this analysis we conclude that 
No source, DETO , HED6 and (marginally) DDTe models are disfavored by the data, while other models are compatible 
with the data. 
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TABLE 8 

Pairwise comparison of models from Table \2\ For each model the maximum 

DEVIATION \pbf/cr p \ WITH RESPECT TO ALL OTHER MODELS IS GIVEN. 


Model 

ISGRI Sc SPI(400-1350 keV) 
Maximum deviation 

ISGRI Sz SPI(70-1350 keV) 
Maximum deviation 

No source 

9.0 

ii.i 

DDTlpl 

0.7 

1.8 

DDTlp4halo 

1.1 

2.1 

DDTe 

2.2 

3.6 

DETQ 

12.5 

12.6 

HED6 

4.3 

5.4 

W7 

0.8 

1.3 

ddtlp4 

2.3 

1.3 

DD4 

1.6 

1.0 


Note. — The ISGRI and SPI data for the early and late period are used jointly. Large 
deviations imply that the model is disfavored by the data. Using a threshold of 2 standard 
deviations one can conclude that No source, DETO , HED6 and (marginally) DDTe models are 
disfavored by the data, while other models are compatible with the data. Since 3Dbbal 
model that was designed to fit the INTEGRAL data in the early period, this model has 
been excluded from this test. 
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